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Abstract 

We have studied the hard magnetic properties of the low-temperature phase of MnBi with 
first principle calculations based on the density functional theory. The calculations have been 
carried out on two distinct unit cell configurations MnBi and BiMn with the element in the 
unit cell origin named first. Our results show that these configurations are not equivalent and 
that MnBi describes the system better near T = OK and the BiMn configuration describes the 
system better for T > 300A7 The magnetic moments of both configurations agree well with 
experimental measurements considering both spin and orbital contributions. At high tempera¬ 
tures the magneto-crystalline anisotropy energy increases with increasing unit cell volume and 
reaches a maximum of 2.3 MJ/m 3 and a c/a ratio of 1.375. 
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1 Introduction 

The low-temperature phase (LTP) of MnBi has interesting properties suitable for permanent 
magnetic applications [T]. At low temperatures its spin configuration is in-plane but at T = 80 K 
a spin reorientation takes place and the magneto-crystalline anisotropy energy (MAE) increases 
with increasing temperature up to 2.1 kJ/m 3 at T = 500AT. At T = 593 K a first-order phase 
transition occurs reducing the c/a ratio [21 - These properties motivate a systematic study 
of the magnetic properties of MnBi in dependence on the change of lattice parameters with 
temperature. 

The LTP phase of MnBi is a hexagonally close packed structure with ABAC stacking 
and described by space group 194 P6 3 /mmc. The A layers have the Wyckoff position 2a the 
alternating B and C layers 2c. The 2 a and 2c positions have different atomic environments, 
which means that exchanging Mn and Bi positions leads to different systems (see fig. [I]). For 
clarity, the configuration with Mn at the 2a site will be called MnBi and the configuration with 
Bi at the 2 a site will be called BiMn throughout the paper. Sources for both configurations 
can be found in literature (e.g. BiMn: [3J|3], MnBi: [3,, 6]) 
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Figure 1: Atomic environments in the MnBi unit cell, (a) Trigonal prism around 2c cite, (b) 
Distorted square anti-prism surrounding 2 a site. Figures created with VESTA [?]. 


2 Method 

Ab-initio calculations based on the density functional theory (DFT) have been performed. 
The central statement of the DFT is that the total energy of a many-electron system can 
be expressed as an functional of the electron density. This leads to the formulation of the 
Kohn-Sham equation, a Schrodinger-like equation on a continuous electron density [3- The 
Kohn-Sham equation can be solved in a self-consistent manner. 

The calculations have been carried out with the WIEN2k package, a full-potential, aug¬ 
mented plane wave code B ESI- As exchange-correlation potential the generalized gradient 
approximation in the formulation of Perdew, Burke and Ernzerhof (PBE-GGA) has been used 
m- The main advantage of GGA over the local (spin) density approximation L(S)DA is that 
not only the electron density, but also its first derivative is taken into account. 

Spin-polarized, scalar-relativistic calculations on both MnBi configurations have been per¬ 
formed until convergence with a criterion of 10 ~ 9 Ry on a 25 x 25 x 15 k-mesli has been reached. 
After that, the relativistic effects such as the orbital contraction of Bi s-states and the spin-orbit 
coupling (SOC) have been introduced. The MAE is the difference of the total energies with 
SOC in easy and hard direction: MAE = Fl[ioo] — F^ooi] ■ It is worth noting that the results 
depend on the SOC initialization method and the resulting crystal symmetries. While the error 
due to incorrect symmetries is small compared to the total energy of the system (10 4 Ry), it 
is in the order of typical MAE values (10 ~ 5 Ry). Orbital contributions to the total magnetic 
moments have been determined afterwards by calculating the density matrix. 


3 Results 

3.1 Lattice Optimization 

Based on the lattice parameters reported in [2j calculations to find the equilibrium lattice 
parameters of MnBi have been performed. Optimizing unit cell volume and c/a ratio yields 
a = 4.3408A and c = 5.7142A with c/a = 1.3164. Both lattice parameters deviate from the 
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experimentally measured values (a = 4.285A, c = 6.113A in 0 ) but are in good agreement 
with other DFT calculations mm 

The equilibrium lattice parameters correspond to T = OAT. The MAE of the MnBi system 
was determined to be —0.11 MJ/m 3 which is in better agreement with experimental mea¬ 
surements (MAE = —0.2 MJ/m 3 at T = 4A", (TJ (2|) compared to other DFT calculations 
(MAE = —2.1 MJ/m 3 , [SHU]). Calculations on the BiMn system wrongfully predict an uniax¬ 
ial anisotropy with MAE = +1.55 MJ/m 3 at T = OK. 

3.2 Magnetic Moments 

The total spin magnetic moment per MnBi unit cell is 7.37 ^b- The main contribution to the 
spin moment come from the Mn atoms (3.77 /.ib) especially from d electrons (3.68p-s). The 
spins of the Bi atom are aligned anti-parallel to the Mn spins lowering the total spin moment. 
The interstitial moments (0.07/1^) are contributions from electrons outside the muffin-tin radius 
and cannot be assigned to a single atom. Detailed per-orbital contributions are summarized in 
table [TJ 

These results are in good agreement with other DFT calculations, but lower than experi¬ 
mental results with jiMn = 3.82 — 4.25^b mm- Considering the orbital moments as well (see 
table [2]) this experimental value can be reached without the introduction of an phenomenologi¬ 
cal Hubbard U potential. The magnetic moments per atom can be increased by 5% by variation 
of the lattice parameters, however, the magnetic moment per volume stays constant. 

Magnetic moments of the BiMn configuration are between 9% and 12% higher than MnBi 
systems throughout the examined range of lattice parameters and agree with the upper bound 
of experimental values (see tables [3] and [4]) . 




s 

P 

d 

f sum 

Mn 


0.05770 

0.02721 

3.68389 

-0.00008 3.76872 

Bi 


0.02100 

-0.15420 

0.01055 

0.00432 -0.11833 

interstitial 




0.07154 

sum 


0.07870 

-0.12699 

3.69444 

0.00424 3.72193 


Table 1: Spin Moments of MnBi at T - 

= 0A' in hb 



P 

d 

f 

sum 


Bi 

-0.00011 

0.07802 

0.00020 

0.07811 


Mn 

-0.00413 

-0.00015 

-0.00010 

-0.00438 


sum 

-0.00424 

0.07787 

0.00010 

0.07373 


Table 2: 

Orbital Moments of MnBi at T = OK in 



s 

P 

d 

f sum 

Bi 


0.02437 

-0.07656 

0.01122 

0.00558 -0.03539 

Mn 


0.08713 

0.04057 

3.81380 

0.00024 3.94174 

interstitial 




0.28596 

sum 


0.11150 

-0.03599 

3.82502 

0.00582 4.19231 


Table 3: Spin Moments of BiMn at T = OK in 
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p 

d 

f 

sum 

Bi 

Mn 

-0.00754 

-0.00171 

-0.00046 

0.10225 

-0.00052 

0.00026 

-0.00852 

0.10080 

sum 

-0.00925 

0.10179 

-0.00026 

0.09228 


Table 4: Orbital Moments of BiMn at T = OK in fig 

3.3 Temperature Dependence of MAE 

The temperature dependence of the magnetic properties have been calculated by variation of the 
lattice parameters. Only thermal expansion has been considered in these calculations, other 
thermal effects like temperature smearing of electron bands or phonons have been ignored. 
Measurements show that in the temperature range from 300/\ to 700 K the c/a ratio varies 
between 1.35 and 1.45 and the unit cell volume increases from 97A' to lOlA . 

The change of lattice parameters is not sufficient to change the sign of MAE of the MnBi 
system, but the values obtained from the BiMn system are in good agreement with experimental 
results (see fig. [2]). With increasing temperature (unit cell volume) the MAE of the BiMn 
system increases. The MAE is maximized with c/a ratios around 1.375 and the influence of the 
c/a ratio is higher at lower volumes. 



c/a 


Figure 2: MAE dependence on unit cell volume and c/a ratio, BiMn configuration. The black 
line shows the temperature dependence of volume, c/a ratio and MAE. In the dotted region 
the first order transition takes place leading to c/a reduction. 


4 Conclusion 

We have successfully calculated the magnetic properties of LTP MnBi with two non-equivalent 
lattice configurations: MnBi (Mn at the origin) describes the properties around OK and BiMn 
(Bi at the origin) describes the high-temperature properties from 300 K to 700K. Both con¬ 
figurations have been confirmed experimentally in literature. The results suggest that ther¬ 
mal expansion alone cannot explain the unusual temperature-dependent magnetic properties of 
MnBi, but they could be a consequence of a transition between MnBi and BiMn configuration. 
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The liigh-temperature MAE increases with increasing unit cell volume and is maximized by c/a 
ratios of 1.375. 
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